ULB-PMIF-92/10 



Bianchi Cosmological Models and 
Gauge Symmetries 



Olivier Coussaert 1 and Marc Henneaux 2 

Faculte des Sciences, Universite Libre de Bruxelles, 
Campus Plaine CP. 231, B-1050 Bruxelles, Belgium 



1 Chercheur IRSIA 

2 Maitre de Recherches au Fonds National de la Recherche Scientifique. Also at 
Centros de Estudios Cientfficos de Santiago, Casilla 16443, Santiago 9, Chile. 



Abstract 



We analyze carefully the problem of gauge symmetries for Bianchi 
models, from both the geometrical and dynamical points of view. 
Some of the geometrical definitions of gauge symmetries (= "homo- 
geneity preserving diffeomorphisms" ) given in the literature do not 
incorporate the crucial feature that local gauge transformations should 
be independent at each point of the manifold of the independent vari- 
ables ( = time for Bianchi models), i.e, should be arbitrarily local- 
izable ( in time). We give a geometrical definition of homogeneity 
preserving diffeomorphisms that does not possess this shortcoming. 
The proposed definition has the futher advantage of coinciding with 
the dynamical definition based on the invariance of the action ( in La- 
grangian or Hamiltonian form) . We explicitly verify the equivalence of 
the Lagrangian covariant phase space with the Hamiltonian reduced 
phase space. Remarks on the use of the Ashtekar variables in Bianchi 
models are also given. 

1 Introduction 

It is a well known fact that the distinction between global and local symme- 
tries is somewhat subtle in the case of spatially homogeoneous cosmological 
models of the Bianchi type [|I], |^, |3| . There are indeed at least three different 
ways to define the gauge symmetries Q: 

(i) geometrically, as those diffeomorphisms that preserve manifest homo- 
geneity; 

(ii) dynamically, from a spacetime covariant point of view, as those trans- 
formations that leave the action invariant and that are local in time; 
and 

(iii) dynamically, from the Hamiltonian point of view, as those canonical 
transformations that are generated by the first class constraints. 

1 We restrict the analysis of this paper to class A Bianchi cosmologies as these are the 
only ones for which one can inject the homogeneous form of the metric in the action prior 
to computing the variationnal derivatives |ij ^[ || 
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It is not a priori clear that these three different points of view should yield 
identical conclusions. Futhermore, while the dynamical approaches based on 
(ii) and (hi) yield unambiguous answers, the geometric approach is not so 
well defined since it depends on what one means exactly by "diffeomorphisms 
preserving manifest homogneity" . 

The purpose of this paper is to show that the three different concepts of 
gauge symmetries are actually the same provided one adopts a natural def- 
inition of "homogeneity preserving diffeomorphisms" . This definition agrees 
with the one given in |], |] but differs from that proposed in ||. We first 
consider the metric formulation of Bianchi models and then briefly comment 
on the formulation in terms of Ashtekar variables. 

2 Geometric definition of gauge symmetries 
2.1 First definition 

Bianchi cosmological models are spacetimes invariant under a three- dimen- 
sional group of isometries G acting simply transitively on three - dimensional 
spacelike hypersurfaces 0. 

In Bianchi models, one is thus given an action of G on the manifold M. , 

G — Diff(M) (1) 
ReG ^h R eDiff(M) 

and the spacetime metric g fulfills 

h* R 9 = 9 (2) 

Since the action of G is part of the given geometric data, it is natural to 
define the subgroup H C Dif f(M) of "homogeneity preserving diffeomor- 
phisms" as the set of diffeomorphisms which preserve the action of G, i.e. 
which commute with each h R ,R G G 

f eH <-> fh R = h R f VR G G (3) 

If one performs a diffeomorphism k that does not belong to H, one can 
detect it by just comparing the action of G before and after k has been 
performed. The group H defined by (^) is the centralizer of G in Dif f(Ai). 
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In order to determine the general solution of (|3]) , we observe that the 
manifold M. is foliated by the three-dimensional hypersurfaces £ of tran- 
sitivity of G. It is convenient to introduce a non-vanishing vector field v 
transverse to these hypersurfaces and such that the one - parameter group 
{expx°v } commutes with G, 

(exp x°v)h R = h R (exp x°v) (4) 

This can be done by picking an arbitrary curve 7 : a; t— > j(x°) everywhere 
transverse to the hypersufaces E. One then constructs the family of curves 
x° 1 — > hftljlxo)). Since G acts simply transitively, there is one and only one 
such curve through each point. One defines v to the tangent vector field 
to the family. The vector field v fulfills (f|) and x° is constant along the 
hypersurfaces of transitivity. 

Let us now choose a basis in the Lie algebra Q of G and take the corre- 
sponding Killing vectors £ a (a = 1,2,3) tangent to the hypersurfaces S. The 
vector fields v and £ a provide a basis of the tangent space at each point of 
the spacetime manifold Ai. One has 



[v,Q=0 (5) 

where C c ah are the structure constants of Q. 

Let X a be a basis of invariant vector fields fulfilling also [f,X a ] = 0, 

[v,X a ] = (6) 

[Za,X b ] =0 

One may take X a = £ a on the curve given above. In that case, on finds in 
addition to §)-(§) 

[X a ,X b ] = -C c ab X c (7) 

The most general "infinitesimal" diffeomorphism commuting with G is 
determined by a vector field £ fulfilling 

[£, 6J = (8) 



2 These are the only ones considered here. 
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If one expands £ in the basis {v, X a }, 

£ = fxv + fx a X a (9) 

one gets 

<9 e /i = 0, <9 e /i 6 = 0, (10) 

i.e., 

£ = fj,(x°)v + fi a (x°)X a (11) 

Thus, a homogeneity preserving diffeomorphism depends on four func- 
tions of x°, irrespectively of the homogeneity group and of the spatial topol- 
ogy. This is quite a satisfactory result, since in the non- homogeneous case, a 
general diffeomorphism depends also on four functions ( which may be non- 
homogeneous). Futhermore, the fact that ( pTT|) involves arbitrary functions of 
x° implies that one can arbitrarily localize the homogeneity preserving dif- 
feomorphisms in time ( there is, of course, no hope to localize them in space 
in a spatially homogeneous context). This localizability property appears to 
be an essential feature of gauge (= "local") transformations. 



We close this subsection by observing that (|Tl|) acts though the adjoint 
action on the invariant frame {v, X a } , 

6X a = %X a ]=-n h C* b X c 



This yields 



r dgoo dg 00 d[i a 



dg 0a dfi dfi d , 

9oa = ^Ix^ + 9oa dx^ 9ab dx^ ^ da9ob ( > 

Sg ab = ^ + ^\C c da g c b + C c db9ac ) 

for the components gx^x ) of the invariant metric in the invariant frame 

{v,x a } 

The transformations (|13|) are irreductible if and only if the adjoint repre- 
sentation is faithful, i.e., if the center of Q vanishes. Otherwise, there exist 
non trivial gauge parameters for which the transformations (O) are zero. 
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These are given by /x = 0, ^ = with /i a a non zero solution of ii a C c ah = 0. 
An irreductible description of the gauge transformations can be obtained as 
follows. Let £4 be a basis of the center of Q (Cf b = 0) and £ a supplementary 
vectors so that {£a)£oJ is a basis of Q. On can replace by the gauge 
parameter e A without violating locality in time ( since fi A does not appear 
in which yields 



cfcc c?x° ' dx° 
x dg 0a dfi B dfi 13 s b 

9 ° a = ^ ~dx^ 9 ° a d^ 9aBE + ^P^+^ C 8a90b (14) 

$9ab = H-^p + n\C c Sa g cb + C c 5b g ac ) 

(For more information on reductible gauge transformations, see for in- 
stance || chapter III). 



2.2 Second definition 

A weaker form of "homogeneity preserving diffeomorphisms" has been pro- 
posed in Pj. These are the diffeomorphisms which map any spatially homo- 
geneous metric on a spatially homogeneous metric. For this to be the case, 
one must have in place of @ 

[£,£.] =J&e (15) 

where the matrix n b a is constant. Clearly, any solutions of @ is also a solution 
of (jl5|). The converse is not true for non compact cosmologies of type J, II , 
VI Q or VIIq || . In finite term, the equation ([15|) becomes 

fhaf' 1 = h R > 

Thus, / is a homogeneity preserving diffeomorphism according to the second 
defintion if and only if its belongs to the normalizer of G in Dif f{M). 

Let us now determine which solutions of (|15|) can be localized in an ar- 
bitrarily small spacetime interval. If £(x°) is zero for x° ^ [x^xSJ], one h- as 
ix b a = since [£, £ a ] = for x° £ [^1,^2] anc ^ * s constant. Hence, the 
only homogeneity preserving diffeomorphisms of || that can differ from the 
identity in a finite time interval are precisely the homogeneity preserving 
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diffeormorphisms of the previous subsection. If one feels that the concept 
of "gauge transformation" ( "of the second kind " as one used to say in the 
old days) should involve in a crucial manner the freedom to change the local 
reference frame independently at each time ( localizability of the transforma- 
tion), then, one must impose the condition = in flT5|) and one falls back 
on the first definition given previously. This point has not be addressed in 



3 Lagrangian definition of gauge symmetries 
3.1 Gauge transformations 

The previous considerations are purely kinematical and do no refer to the 
action principle. To decide the issue of which geometrical definition is prefer- 
able, one should introduce physical considerations, namely, study the trans- 
formation properties of the action principle. 

The equations for the Bianchi models of class A can be derived from the 
Hilbert action in which the symmetry of the metric is enforced before one 
computes the variationnal derivatives ( "reduced variationnal principle" ) || 



The lagrangian contains reference to the homogeneity group through the 
structure constants C% c that appear as parameters. It explicit form is[|, [| 



I- 



S[g ab ,N(x°),N a (x°)} = Jdx L( 9l 



dg ab 
dx° 



N, N a ) 



(16) 



L(9ab, 



dg a b 
dx° 



, N, N a ) = J dx 3 N^(R - K c c K d d + K d c K c d ) (17) 



where: 



7~> foS~ic r^o- d i /^idca/^i i r~icdar~i 

ti — — — cd^ a + t> <^cda T ^cda 



K ab = ^i(C ab c + C ba c)N c -^} 



with: 




G 



abc 



bd cerya 
9 9 ^ de 
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Note that the shift iV a appears only through the combinaision N a C^ b , so that 
the components of N a in the center drop out. 

A gauge transformation of ( |TED is a transformation of the dynamical vari- 
ables which must fulfills two conditions ||: 

(i) it leaves the action invariant ( up to a boundary term); and 

(ii) it is localizable; i.e., it must depend on arbitrary functions ofthe coor- 
dinates, which can be chosen so that the transformations reduce to the 
identity outside a finite region. In the present case where we impose 
spatial homogeneity, the independent coordinates reduce to x° and a 
finite region is typically a finite interval [ 

It is the second point that makes the distinction between gauge ("local") 
transformation and global symmetry transformations. Gauge symmetries 
lead to non uniqueness of the Cauchy problem while global symmetries do 
not. 

It should be stressed that the definition (i)-(ii) of the gauge symmetries 
isentirely "covariant" i.e., formulated from a spacetime point of view. The 
fact that x° plays a priviledged role is a result of the spatial homogeneity 
alone and not not of the Hamiltonian splitting of spacetime into space and 
time. 

It is easy to verify that the transformations (|T3| ) or ( |14D when rewritten 
in terms of g a b , the lapse and the shift, leave the action invariant. 
Furthermore, as we have discussed, they are local in the sense of (ii). Thus, 
they define gauge symmetries. They actually provide a complete set of gauge 
symmetries, because once the gauge freedom (^) is frozen- say by imposingg 
N = l,N a = 0- the equations of motion following from (|16|) have a unique 
solution. There is no futher degeneracy besides the one implied by (|TTD. [For 
more on the general description of gauge symmetries, see || ]. Accordingly, 
there is complete agreement between the first geometric definition of gauge 
symmetries and the dynamical definition based on the Lagrangian. 

The extra symmetries which the second geometrical definition allows (el- 
ements of the normalizer of G which are not in the centralizer) fail to be 
dynamical gauge symmetries because, as we have already emphasized, they 
do not fulfill (ii). Hence, they lead neither to a degeneracy of the Lagrangian 
nor to an indeterminacy of the Cauchy problem. Trying to "gauge fix" them 
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is unnecessary to get an unambiguous evolution and would remove true de- 
grees of freedom. Even though these symmetries are not gauge symmetries, 
they do define, however, global symmetries as we now discuss. 



3.2 Global symmetries 



The transformations Q15D with n° a ^ induce the following changes of the 
dynamical variables , 

5N = 

5N a = a a h N b (18) 

$9ab = —O'bdcb — O'bdac 

where a£ is a constant matrix belonging to the Lie algebra of the automor- 



phism group of Q. The transformations ( jig) leave the action (|T6|) invariant 
provided Tr a = and hence, define global symmetries when Tr a = 0. By 
Noether theorem, they give non vanishing conserved charges. 

A crucial feature of the transformations ([Tj|) of the dynamical variables 
is that they are well defined even if the spacial sections are compact. In that 
case, they cannot be induced by a globally defined diffeomorphism[|[, but 
still, they make perfect sense. Futhermore, they take exactly the same form 
irrespectively of the topology. From the dynamical point of view, they should 
always be regarded as global symmetries since they cannot be localized in 
x° by a choice of the parameters a%. In the type I case, the corresponding 
conserved charges account for the integrability of the model |], |2|. 

3.3 Counting degrees of freedom- Covariant phase 
space 

The number of degrees of freedom of a gauge theory can be counted directly 
from the Lagrangian, without having to go through the Hamiltonian analysis. 
It is equal to (half) the dimension of the covariant phase space, defined 
to be the space of solutions of the equations of motion modulo the gauge 
transformations || [TO], [OJ p| . 



Now, a solution of the equation of motion is completely characterized by: 
(i) a choice of lapse N(x°) and shift iV a (:r ); and (ii) initial data g a b{%° — 0), 
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^jztr(x° = 0) at x° = (say) fulfilling the constraint equation 

nn J(N 2 - N S N S ) 
G = <R + K\K] - K)K{) (19) 

2 ■ N i 

G 0i = jf{-C l iK' s + C\ 3 K^) + —(R + KIK]- KyK**) 

The number of independent equations ( P^| ) is equal to 4 — dim{centerQ). 
By a choice of gauge, on can take N — 1 and N a = 0. The residual gauge 
transformations acting on the initial data g^, -4% take the form 

Sg* = ^ + li s (Csa9* + C c Sb g ac )(20) 
x dgab r „„ 1 dg cd d dg ab dg u kl dg ak s dg cb dg a c 

and involve effectively 4—dim(centerQ) parameters ( recall that \i a = (/i A , /i 5 ) 
where \i A is in the center). Hence, the dimension of the space of solutions 
of the equations of motion modulo the gauge-transformation is equal to 12 
(number of g ab and N and iV a do not count since they can be "gauged 
away") minus 4 - dim(center Q) ( number of independent transformations 
(120]) . The number of degreees of freedom is half this count, 

n = ^{12 - 2 [4 - dim(center Q)]} (21) 
= 2 + dim(center Q) 

The count fl2jj ) agrees with the count of 0, Chapter 3, where it is shown 
that n is equal to the number of dynamical variables ( here 10 = number of 
g^y) minus the number of times the gauge transformations strike. A gauge 
transformation strikes n times if the corresponding gauge parameter together 
with its time derivatives up to order n-1 appear effectively in the gauge 
transformations]. Here, \i and /i a strike twice , but the gauge parameters /i A 
in the center strike only once. Hence, again, 

n = 10 — 2{1 + 3 — [dim(center £/)]} — dim(center Q) (22) 
= 2 + dim(center Q) (23) 
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The symplectic structure on the covariant phase space is determined by 
the boundary term in the variation in the action. One finds here (on-shell) 

5S = [U ab 5g ab \ll (24) 

with 

Tl ab = Jg(g ab K\ - K ab ) (25) 

Thus, the symplectic two-form on the covariant phase space is the two- form 
induced by the canonical two-form 

m ab A 8g ab (26) 

on the quotient of the surface ( |T9| ) by the transformations (|20|). [ One easily 
verifies that the induced two-form on the surface (|19|) has a rank equal to 
4 — dim(centerQ) , with zero eigenvectors determined by (|20|). It is thus 
regular on the quotient space.] 



4 Hamiltonian analysis 

The reader familiar with the hamiltonian analysis of Bianchi models will 
certainly have noticed that the above results, obtained from a spacetime 
covariant point of view, are just the same as the Hamiltonian results (see 
e.g. 0, |!|). The U ab are the momenta conjugate to g a b, the constraints 
( |19|) are the Hamiltonian constraints and are first class. The transformations 
( |2"0"D are just the transformations which they generate. The details may be 
found in the literature and will not repeated here. We can thus conclude 
not only that the various concepts of gauge transformations do coincide, but 
also, that the covariant phase space is identical (as symplectic manifold) with 
the Hamiltonian reduced phase space. This later point is of course not an 



accident since it is been established on general grounds in [O] 



5 Coset spaces 

One can generalize Bianchi models by assuming local homogeneity instead 
of global homogeneity. In the globally homogeneous case, the surface of 
transitivity can be identified with the homogeneity group G. In the locally 
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homogeneous case, one replaces the surfaces of transitivity G by the quotient 
space G/K where K is a discrete subgroup of G acting on G properly dis- 
continously. The diffeomorphisms of G that do not commute with K fail to 
globally defined on the quotient space. Homogeneous objects do remain well- 
defined, however. To study coset spaces, it is convenient to lift everything 
from G/K to G. One then applies the techniques of the previous sections. 
One finally takes the quotient at the very end. Since the gauge transfor- 
mations are generated by the vector fields X a , which commute with £ a and 
hence with G and K C G, they are well defined in the quotient. 



6 Remarks on Bianchi models and Ashtekar 
variables 

Recently, Ashtekar has introduced new variables in terms of which the con- 
straints of general relativity take a polynomial form ( see e. g. the reviews 
[14,15]). In order to gain insight into those variables, various authors have 
reformulated Bianchi models in terms of them [|l^, [17] . An interesting ques- 



tion then arises: what is the most general form of the Ashtekar variables 
compatible with homogeneity? We first address this question. 

A SO(3)-connection A^ k and a triad e/ a )k are homogeneous if and only 
if the images h* R A^ a)k and h* R e^ a)k of A^ a)k and e (a)k differ from A^ a)k and e (a)fe 
by a local SO (3)- transformation JT8|, |T§, |2D|, 0} ^2 



h* R A k = i)A k ^) R l + i) R d k i) R 
h* R e k = ip R e k 



(27) 
(28) 



where ip R is a spacerime-dependent element of SO (3) that depend also on 
R G G . If we assume that the triads are non-degenerate, or degenerate at 



most at isolated points, one infers from 



81) and h 



R1R2 



h* Rl h* R2 that 



(29) 



so that the difficulties pointed out in |23[ and characteristic of Abelian gauge 
groups do not occur. One can then apply the general techniques developed 
23], 2~5| . Since the group G acts simply transitively, one may define x{P) 



111 

as 



X (P)=^ R \P ) 



(30) 
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where Pq is a fiducial reference point on the surfaces of transitivity to which 
P belongs ( for instance, one can take it to be the intersection point of S 
with the curve 7 ) and R is the unique element of G such that h R (P ) = P. 

Let us define A and e by making the S0(3) gauge transformation generated 
by X, 

A = xAx' 1 + Xdx' 1 (31) 
e = xe 

There, one finds from (p7|) , (p8|) , (p9|) and (^) that A and e are strictly in- 
variant, 

h* R A = A (32) 
h* R e = e (33) 

There is thus no loss of generality in assuming the Ashtekar variables to 
be strictly homogrneous. This confirms by a geometrical method the findings 
of |n| |17| . In the invariant frames {A a }, the components of Ar a ^ and e(a)k 
depend only on x° . The Einstein equations define a problem of classical 
mechanics, as in the metric formulation. 

The definition of the gauge symmetries proceeds then as above. These 
symmetries are found to include besides the homogeneity preserving diffeo- 
morphisms of the metric formulation , the local SO (3) internal transfor- 
mations acting on the internal indices. There is again complete agreement 
between the first geometrical definition and the dynamical definitions. 

The analysis can be extended to the discussion of discrete symmetries (see 
2(| |27||) provided we assume the connection to be an 0(3)-connection rather 



that an SO (3) -connect ion . To illustrate the point, let us consider Bianchi 
type I models, ds 2 = —N 2 (x°)(dx ) 2 + g a b(x°)dx a dx b , with the discrete sym- 
metry a; 3 1 — > —x 3 . Then g 13 = g 2 3 = 0. In terms of the triad e^k{x°) and 
the connection co>( a )fc(x°), the condition of invariance becomes 

e k ia) = X (a \ b) ei b) (34) 

-r=x {a \ b ^+^ iabc) x(M ) d k x~ i u (35) 

where e' fe (a) = (ef^e^, -e { 3 a) ) and u'jf = (wJ o) ,^ o) , -W3 ) 
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Equation ( p^) shows that X n>) ^ as ^ wo eigenvalues equal to +1 and one 
equal to -1. Hence it does not belong to SO (3) but rather to 0(3). In this 
case, one must extend ( |3~3| ) in the following way: 

4 (a) =X (a) (6) 4 6) (36) 
ui'^-X^P + leWxwdkX- 1 ® (37) 

for elements X n>) °f 0(3) which are not in the connected component to the 
identity. Indeed, u'ff^ is a pseudo- vector for 0(3) transformations. 
In an appropriate basis, the matrix X n>) rea ds 

(1 \ 
X {a \ b) =010 (38) 

\0 -i I 

and is constant. The invariance condition implies then 

e« = e< 2) = 0, eS 3) = ef = (39) 
u m 1= JV = , UJ ? ) = 4 2) = (40) 

The variables remaining in the reduced variationnal principle for Bianchi 
type 1_ cosmologies with discrete symmetry x 3 i— > — x 3 are thus: 

(i) ei a) ,4 3) (a = 1,2; A; = 1,2), 

(ii) ef,4 fc) (« = 1,2), 

(iii) iV and iV fc (which drops out). 

The gauge symmetries are (i) the redefinition of the time variable; and (ii) 
the SO(2)-rotations in the (l)-(2) plane. 



7 Conclusion 

We have shown that there exists a natural definition of homogeneity pre- 
serving diffeomorphisms for Bianchi models such that none of the difficulties 
found in [3] arises. With that definition, the geometrical and dynamical con- 
cepts of gauge symmetries coincide irrespectively of the spatial topology. We 
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have also verified the isomorphism of the Lagrangian covariant phase space 
with the Hamiltonian reduced phase space. Finally, we have discussed spa- 
tially homogeneous Ashtekar variables. We have illustrated how the general 
theory of invariant Yang-Mills connections apply. This theory is actually 
relevant in the discussion of all minisuperspace models formulated in terms 
of Ashtekar variables, and not just the Bianchi cosmological models with a 
three-dimensional simply transitive invariance group. In particular, one can 
generalize the considerations of [28] section V to the Ashtekar description. 



8 Acknowledgements 

We are grateful to Luca Bombelli for useful discussions. This work has been 
supported an part by the FNRS (Belgium) and the Commision of the Euro- 
pean Communities. 

References 

[1] M.Henneaux, "Sur la geometrodynamique avec les champs spinoriels", 
chapter V, PhD Thesis (Brussels 1980). 

[2] M.Henneaux, Ann. Inst. Henri Poincare 34(1981)329; see also 
M.Henneaux, Lett. Nuovo Cimento 38 (1983) 609 (especially page 612). 

[3] A. Ashtekar and J.Samuel, Class. Quant. Grav.8(1991)2191. 

[4] R.T.Jantzen, Commun.Math.Phys.64(1979)211. 

[5] M. A. H.Mac Callum, in "General Relativity: An Einstein Centenary 
Survey", S.W. Hawking and W. Israel eds, Cambridge University Press 
(Cambridge: 1979). 

[6] G.E.Sneddon, J.Phys.A9(1975)229. 

[7] A.H.Taub, Ann.Math.53(1951)472. 

[8] M.Henneaux and C.Teitelboim, "Quantization of Gauge Systems", 
Princeton University Press(Princeton: 1992). 



14 



[9] J.Kijowski and W.Szczyrba, Comm.Math.Phys.46 (1976) 183. 
[10] E.Witten, Nucl.Phys. B276 (1986)291. 
[11] G.Zuckerman, Yale University Preprint(1986). 

[12] A.Ashtekar, L.Bombelli and O.Reula, in "Analysis, Geometry and Me- 
chanics: 200 years after Lagrange", eds.M.Francaviglia and D.Holm, 
Noth Holland, Amsterdam 1990. 

[13] G.Barnich, M.Henneaux and C.Schomblond, Phys.Rev.D44(1991)R939. 

[14] A.Ashtekar, Lectures on Non-Perturbative Canonical Gravity, World 
Scientific, Singapore (1991). 

[15] C.Rovelli, Class.Quant.Grav.8(1991)1613. 

[16] H.Kodama, Prog.Theor.Phys., 80(1988)1024. 

[17] A.Ashtekar and J.Pullin, Ami.Israel.Phys.Soc. 9(1990)65 

[18] PCordero and C.Teitelboim, Ann.Phys.(N.Y.) 100(1976)603; 
R.Benguria,P.Cordero and C.Teitelboim, Nucl.Phys.B122(1977)61. 

[19] V.N.Romano,A.S.Schwarz and Yn.S.Tyupkin, Nucl.Phys. 
B130(1977)209. 

[20] P.G.Bergmann and E.J.Flaherty, J.Math.Physl9(1978)212. 

[21] A.Trautman, J.Phvs.A13(1980)Ll. 

[22] A.S.Schwarz, Comm.Math.Phys.56(1977)79. 

[23] M.Henneaux, J.Math.Phys.23(1982)830. 

[24] J.Harnad,S.Shnider and L.Vinet, J.Math.Phys.20(1979)931. 

[25] PForgacs and N.S. Manton, Comm.Math.Phys.56(1977)79. 

[26] B.G.Schmidt, Comm.Math.Phys.l5(1970)329. 

[27] M. A. H.Mac Callum, in Cargese lectures in Physics volume 6, 
ed.E.Schatzman, Gordan and Beach (New York 1973). 

15 



[28] K.Kuchar, Phys.Rev.D4( 1971) 955. 



16 



